ON A PROBLEM OF CHEN AND LIU CONCERNING THE 
PRIME POWER FACTORIZATION OF 71! 
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Abstract. For a fixed prime p, let e p (n!) denote the order of p in the prime 
factorization of n\. Chen and Liu (2007) asked whether for any fixed m, one 
has {e p (rt 2 !) mod m : n G Z} = Z m and {e p (q\) mod m : q prime} = Z ra . We 
answer these two questions and show asymptotic formulas for #{n < x : n = 
a mod d, e p (n 2 \) = r mod m} and #{q < x : q prime, q = a mod d, e p (q\) = 
r mod m}. Furthermore, we show that for each h > 3, we have {n < x : n = 
a mod d, e p (n h \) = r mod m} 3; 4 /( 3 ' 1 + 1 ). 



1. Introduction 

Let pi = 2, P2 = 3, ... be the sequence of prime numbers in ascending order and 
consider the prime factorization of 

ne p (n!) 

Pj ^ n 

Legendre $M p. 10-12] (see also p. 263], [SI Ch. 1.3]) showed that for any 
nonncgative integer n and any fixed prime p we have 

(1.1) e p (n!)=]T 

i>i 

where s p (n) denotes the sum of the digits of n in base p, i.e., 

where £i(n) € {0, 1, . . . ,p — 1}. A well-known area of application for e p (n\) is the 
determination of the explicit numerical error term in Mertens first theorem |161 
Ch. 1.4]. The investigation of the distribution properties of e p (n\) can be said to 
have started with Erdos and Graham |5J p. 77] who stated (in our notation) that 
"it is annoying that we cannot even show that for all k there is an so that in 
the prime decomposition of n^l all the e Pj (rik^) 7 1 J ' *S k, are even." In 1997, 
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Berend [I] solved this problem by showing that for any fixed m > 2 there are 
infinitely many n that satisfy 

e Pl (n!) = e P2 (n!) = • • • = e Pk (nl) = mod ra, 

and the set of all such n has bounded gaps. In his solution, Berend [1] strengthened 
the problem of Erdos and Graham in two different directions. On the one hand, he 
not only considered the parity of the exponents but studied more generally if they 
were divisible by a fixed integer m ^ 2. On the other hand, he already treated 
subsets of integers with prescribed multiplicative properties instead of looking at 
the entire set of integers n. In particular, he showed that for arbitrary fixed positive 
D, k and m there exist infinitely many n such that all the exponents e Pj ((dn)!), 1 ^ 
j ^ k, 1 ^ d ^ D, are divisible by m. 

Several authors considered in the last years extensions of the Erdos-Graham 
problem, namely, Berend/Kolesnik [2], Chen [3], Chen/Liu [U [5], Chen/Zhu [6], 
Luca/Stanica [IT], Sander [14] and Zhai [17] • The most general result is due to 
Berend and Kolesnik [2] who proved unconditionally that 

^{0 ^ n < x : ii = a mod d, e qj (n\) = Tj mod rrij, 1 ^ j ' ^ k} 



x 



O (x 1 ^ 



dmxmi ■ ■ ■ rrik 

for any integer a and d ^ 1 where k > 1 is fixed, q = (<?i, <72, ■ • ■ , <7fc) is a vector 
of distinct, not necessarily ordered primes, m = (mi, 7712, • • • , rrik) is a vector of 
arbitrary integers > 2, and r = (ri, r%, . . . , rk) is such that Q ^ Vj < rrij for 
j = 1, 2, . . . , k, and <5 = 5(m, q, r) > is effectively computable. 

Intriguing problems arise when the sequence of integers n lying in a fixed residue 
class is replaced by sparser sequences such as primes, squares or higher-degree 
powers. Chen and Liu [5] posed several problems in that respect (see also [17] for 
generalizations of these problems). In particular, at the end of their paper they 
remark that they even have no answer to the following basic questions: 

Question 1: Is it true that for all fixed p and m, 

{e p (n 2 !) mod m : n 6 Z} = Z m ? 
Question 2: Is it true that for all fixed p and m, 

{e p (ql) mod m : q prime} = Z m ? 

Zhai 17\ Theorems 3 and 4] obtained a partial answer to Question 1. He showed 
that for all h > 2 and r 6 Z, there are infinitely many n such that e p (n' 1 !) = 
r mod m provided that 



(1-2) P> 



4m - 2, if h = 2, 
h h m h -\ if/i>3. 



From his proof one can obtain a lower bound of the forrrQ 

#{rt < a; : e p (n h \) = r mod m} ^S> p ,h,m log 2:, £ — »■ 00. 

Unfortunately, Zhai's method cannot be applied in the case of small p, such as to 
treat e2(n 2 l) or e^(n 2 \). 



^By / <Ccj g rcsp. / g we mean that there exists a constant C depending at most on ui 
such that / ^ Cg resp. / > Cg. 
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The aim of the present paper is to use our current knowledge of the distribution 
properties of the sum-of-digits function to give complete answers to Questions 1 
and 2. We are able to improve on Zhai's result and to generalize Chen and Liu's 
questions in two different respects. First, we are able to drop the superfluous 
condition (|1.2p and to find asymptotic formulas for the counting functions in the 
case of squares and primes. Second, we give a general lower bound for h > 3. 

Using our results we get the following nice application: Let Z(n) be the number 
of ending O's in base 10 of n\. Observe that 

Z(n) = min{e 2 (n!), e 5 (n!)} = e 5 (n!). 

Then it will follow from Theorem 12.11 that 

lim — #{n < x : Z(n 2 ) = r mod to} = — 
x— toe x to 

for every to ^ 2 and ^ r < to. The analogous result holds also true for the 

number of ending O's of factorials of primes. 

2. Main results 

In the sequel, let tt(x; a, d) be the number of primes = a mod d that are less than 
or equal to x. 

Theorem 2.1. Let p be a prime, m, d > 1 and a < d, ^ r < to. Then there 
exist constants 5 p ]m > and 8 p 2 ) n > such that 

#{n < x : n = a mod d, e p (n 2 !) = r mod to} 

= £ +0 ((i„ g x)"/ V -a), 

and 

H={q < x : q prime, q = a mod d, e p (q\) = r mod to} 

= V ' ; + O [ (logx V^.™ ) . 

TO V / 

The implied constants depend only on p. 

The proof of this result is notably based on recent work by Mauduit and Rivat [T3] 
and Martin, Mauduit and Rivat [12 , and uses exponential sum estimates of hybrid 
type. In contrast, we use an idea of Stoll [15] to obtain general lower bounds for 
higher-degree powers. The method is constructive. 

Theorem 2.2. Let h > 2, p be a prime, m, d > 1 and ^ a < d, ^ r < to. 

Then, as x — > oo, 

(2.1) #{n < x :n = a mod d, e p (n h l) = r mod to} > p ,h irf , m x 4/(3,l+:L) . 
Moreover, there is an effectively computable constant C = C(p,h,d,m) such that 

{e p (n h \) mod to: 0^n<C, n = a mod d} = 1 m . 

The constant C can be directly obtained from the proof. We remark that 

(2.2) {e p {n h l) mod to : < n < p 1/h + (to - 2)d, n = a mod d} ^ Z m . 

By a probabilistic argument one might expect that we have the full set of residues 
after about to log to steps. However, as (|2.2j) shows, this is not true since there is a 
crucial dependency of p in the bound for n in (|2.2[) . 
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3. Proof of Theorem 12.11 

Legendre's formula shows that 

(3.1) e p (n!) = r mod m n — s p (n) = r(p — 1) mod (p — l)m. 

In order to prove Theorem 12. 1[ we need some auxiliary results. In particular, we 
have to deal with exponential sums containing the sum-of-digits function of primes 
and squares. Let us(b) denote the number of different prime divisors of b. The first 
proposition is a generalization of |13[ Theorem 1] and the second proposition is 
taken from [T21 Proposition 4] . 

Proposition 3.1. Let 6^2 and a, /3, 7 real numbers such that (b — l)a ^ Z. Then 
there exists a constant o^ 1 2 > such that^ 



V" e(as b (n 2 ) + (3n 2 + ~/n) < b (logx) 



(u(b) + lQ)/4 x l-^> 



Proof. This result can be proven in the same way as [T3j Theorem 1]. Thus, we 
just give a short outline. Let b"^ 1 < x < b v and set f(n) = asb(n). As in the 
Mauduit-Rivat case, it suffices to show that 



(w(6)+6)/4 a .l- CT 



(3.2) S x := J2 e (/( n2 ) + ^ + ^ <<b ( lo S x ) 

b u ~ l -Cn^x 

for some constant o~^ a . Lemma 15 from [13j (a van der Corput-type inequality) 
implies that Si is bounded by (some constant times) 



b u- P /2 + b u/2 max 

l^|r|<fcP 



X] e (f(( n + r ) 2 ) +(3(n + r) 2 +7(n + r)) 

1/2 

c(-f(n 2 ) - (3n 2 - jn) 



<h b u - p ' 2 + b^ 2 max 

l^|r|<hP 



^ e(/((n + r) 2 )- f{n 2 ) + 2finr) 



1/2 



where l^p^i//2isan integer which we will choose later on. Set A := v + 2p + 1. 
Using [13l Lemma 16], we obtain 



(3.3) 
where 



Si « h V-p' 2 + b^ 2 max |S 2 | 1/2 

l^|r|<6P 



e(h((n + r) 2 )-h(n 2 ) + 2(5nr), 



and f\(n) is defined by 



f x (n) := a ^ £j(n), 



2 If /3 



0, 1 131 Theorem 1] shows this result with an error term of the form 



(log x)^^ + ^/ 2 x x CT instead of (log x)^^ 1 ^ 10 '/ 4 ^ 1 CT . However, we want to remark that the 
proof given in |13| already implies the better error term as stated in this proposition. 
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where Sj(n) denotes the j-th digit of n. Note, that fx(n) (a so-called truncated 
sum of digits function) sums up just the A lower placed digits (multiplied with a). 
Lemma 17 from |13j (again a van der Corput-type inequality) implies now that 

(3.4) \S 2 \ 2 <b 2 "- 2 ? + b» max |S 3 |, 

l^|s|<& 2 c 

where 

■e(-f x ((n + r) 2 ) + f x (n 2 )-2(3nr) 1 

the interval I(v, s, p) is given by I(v, s, p) = {n G N : b^^ 1 < n,n + sb p ^ b"} and 
p is an integer satisfying 1 ^ p ^ u — 2p — 1. Thus we get that IS3I is equal to 

Note, that the terms containing /3 and 7 are vanished. Mauduit and Rivat consid- 
ered exactly the term £3 and they showed that 

(3.5) |£g| <f, V "(b)+6 b u-2p 

for every 1 < p < v/2, 1 < p sC v - 2p - 1, 1 s$ \r\ < b p and 1 ^ \s\ < b 2p (see [H 
Eq. (45)]). Equations (O, (]33), and dUS]) finally imply 

£1 « fc vWM/iy-p/*, 

As in [13] . it is now possible to choose p and p in order to obtain (|3.2[) . This finishes 
the proof of Proposition 13. II □ 

Proposition 3.2. Lei b ^ 2 and a, /? reed numbers such that (b — l)a Z. TTien 

(2) 

£/iere exists a constant o~, ' > suc/i £/iai 



^ e{as b {q) +pq) < fc (logx)V~ 



<2<a: 

5 prime 



Proof of Theorem \2.1\ We just give a proof of the stated result for the squares. The 
case e p (q\), q prime, can be shown exactly the same way but using Proposition 13.21 
instead of Proposition 13.11 In the following we use the abbreviation 



to' = (p — 1)to. 

Relation (13.11) allows us to write 



< x : n = a mod d, e p (n 2 l) = r mod m} = Tj(x), 

where 

Tj(x) := #{n < x : n = a mod d, n 2 = j mod to', s p (n 2 ) = j — r(p — 1) mod to'}. 
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Using discrete Fourier analysis, we have 



r i(«) = ElJ E 



n<x \ 0^u<d 



e U 



1 \ ( , n 2 — j 



(to/ E e ( ^ 



_L y e fc- - 

to' ' \ ml 

0^fc<m' x 

«p(^ 2 ) - (j - - 1)) 



This can be written as 
1 



n<x 0^u<d 



w-sb-E E M-V 1 ) E -l*^ 



9 

n — J 



0<€i<p-l 



5] ef&m + fc) 



s P (n 2 ) - (j - r{p - 1)) 



and we obtain (splitting the part coming from £2 = and £2 > 0) 



?H*) = i E E E 



n=a mod 



fcn 2 — kj — £\mj + £ims p (n 2 ) 



O 



E E 



dm' 2 



\ 



0^k<m' Q<£ 2 <m 



E 



s p (n ) H -n + -n 

ml m a 



Thus we get that #{w < x : n = a mod d, e p (n 2 !) = r mod m} is given by 
/ \ 



(3.6) MT + O 



where 



max 

0^fe<m' 
. 0<fi<p-l 
\ 0<£ 2 <m 



fcn 2 — kj — £\mj + £ims p {n 2 ) 



0^j<m' n<x 0^k<m' 0^£i<p-l 
ri— a mod d 



Next we calculate the main term MT in (|3 . 6ft . Therefore, let us define lj(n) for 
all < j" < m' and for all positive integer n by 



1, if n = j mod to', 
0, otherwise. 

Then we get that the main term MT is equal to 

1 v ^-^ ^-^ / — £±mj + £ims p (n 2 ) \ 1 



7 E E E 

n—a mod d 



7 E E E e(* 

n=a mod d 



- E ° ^ 



0^fc<m' 



s p(" 2 ) - J 
p-1 



l,-(n 2 ). 
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Since s p (n 2 ) = j mod (p — 1) if tt, 2 = j mod m', we obtain that the remaining 
exponential part is equal to 1 for all nonzero summands. Furthermore, the relation 

E 1 i(n) = l 
holds trivially for any integer n. Thus we finally have 

MT -h E E Un 2 ) 

OsS^i<p-l _n<x 0^j<m' 



n— a mo 



d d 



" 1 E E m-V^+om 



n<a; 0^j'<m' 
n=a mod d 



It remains to bound the error term in (|3.6j) . Since < ^ < m, we have (4m + 
l<x)lm! ■ (p — 1) = i\ + £2/ 'tti ^ Z. Thus we can employ Proposition 13. II (note, that 
u)(p) = 1 since p is prime). Setting 

:= min o^}, , „ \ 1 
p ' 0^i<p-l P.(*i"»+^)/m" 

0<^2<m 

we finally obtain the desired result. □ 



4. Proof of Theorem 12.21 

Proof of Theorem \2.2{ Consider the polynomial t[x) € Z[x] with 

t{x) = d(p — l)m • (m^x 3 + m2i 2 — mix + too) + a, 

where 7713, 777,2, mi, too are positive integers. Lemma 2.1 in |15) says that for all 
u > 1 and 

< mo + -jp ° < p u , 
u(p — ljm 

p u_1 ^ m 2 ,m 3 <p u , 

(4.1) 1 < mi < p u l (hp(6p) h ) 

the polynomial (t(x)) h = Xa=o Ci2; ' e ^"M nas a ^ positive integral coefficients 
with the only exception of the coefficient of x 1 which is negative. Also, note that 
a ^ d(p — 1)to and thus we have 

\d\ < (Ap u d(p - l)m) h . 

In order to have that the range (|4.ip for mi contains at least one admissible integer 
mi we suppose now that u is such that 

(4.2) p u > hp(6p) h . 
Furthermore, let k be such that 

(4.3) p k > (4p u d(p - l)m) h . 

Note that p k ^> p ,h,d, m P uh as u — > 00. We get as in [TS] that 

Sp ((t(p k )) h ) = k(p-1) + M, 
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where M does not depend on k provided k is such as in (|4.3|) . In addition, we have 
t(p k ) = a mod d and {t{p k )) h = a h mod (p — l)m. Therefore, by (|3.1|) . for each k 
with (|4.3p and j > we have 

e pm ^)n = (*(p M - j )) h -^(p fc+J )) h ) 

«' 1 M 

— (fc + j) mod m. 



p — 1 p — 1 

Note that (p — 1) | (a h — M) so that for each fixed r with ^ r < m there is j with 
^ j < m — 1 such that 

e p ((t(p k+ i)) h !) =r modm. 

By construction we thus find ^ p ,h,d,m P Au distinct integers that are all <^ p ,h,d,m 
pu(3h+i) ^ or more details we refer to [H]). This proves (12. ip . 

To get an explicit bound for C (p, h, d, m) we only have to make some admissible 
choices, say, uq and ko, for u in (|4.2j) resp. fc in (|4.3|) . and estimate t(p k+m ^ 1 ). 
First we take uo to be such that p u ° ^ hp 2 (6p) h . Secondly, we find ko such that 
pfco <^ (4:p u °d(p — l)m) h p. It is now a straightforward calculation to find 

C ^ t{p k+m - x ) s: 2(p - l)mp Uo p 3{ko+m -V 

< m 3h+1 d 3h (p - l)^+y™ (Ahp 2 (6p) h ) 1+3h . 
This concludes the proof of Theorem 12.21 □ 

5. Concluding remarks 

We end our discussion with a few remarks. It seems possible to use the approach 
of Drmota, Mauduit and Rivat [5] to get an asymptotic formula in Theorem 12.21 
provided that p is a very large prime whose size is about exponential in h. For 
small p it is already an open and surely very difficult problem to find an asymptotic 
formula for e p (n 3 ) in arithmetic progressions. As a further remark, we also stress the 
fact that Theorem 12 . 1 1 and Theorem 12 .21 hold for arbitrary quadratic polynomials in 
place of n 2 , respectively, for arbitrary P(x) £ Z[x] of degree h (with P(N) C N) in 
place of x h . A minor variation of the used arguments will yield these generalizations. 

References 

[1] D. Berend, On the parity of exponents in the factorization of n\, J. Number Theory 64 (1997), 
no. 1, 13-19. 

[2] D. Berend, G. Kolesnik, Regularity of patterns in the factorization of n!, J. Number Theory 

124 (2007), no. 1, 181-192. 
[3] Y.-G. Chen, On the parity of exponents in the standard factorization of nl, J. Number Theory 

100 (2003), no. 2, 326-331. 
[4] Y.-G. Chen, W. Liu, On the exponents modulo 3 in the standard factorisation of n\, Bull. 

Austral. Math. Soc. 73 (2006), no. 3, 329-334. 
[5] Y.-G. Chen, W. Liu, On the prime power factorization of n!, II, J. Number Theory 122 (2007), 

no. 2, 290-300. 

[6] Y.-G. Chen, Y.-C. Zhu, On the prime power factorization of n!, J. Number Theory 82 (2000), 
no. 1, 1-11. 

[7] L. E. Dickson, History of the theory of numbers, Vol. I: Divisibility and primality, Carnegie 
Inst, of Washington, Washington, D.C., 1919; republication by Dover, Mineola, 2005. 

[8] M. Drmota, C. Mauduit, J. Rivat, The sum-of-digits function of polynomial sequences, J. 
London Math. Soc. 84 (2011), 81-102. 



ON A PROBLEM OF CHEN AND LIU 



9 



[9] P. Erdos, R. L. Graham, Old and new problems and results in combinatorial number the- 
ory, Monographics dc L'Enscigncmcnt Mathematique, 28. Universite dc Geneve, L'Enscigncmcnt 
Mathematique, Geneva, 1980. 

[10] A.-M. Legendre, Theorie des nombres, 3e edition, tome 1, Firmin Didot Frcres, Paris, 1830. 

[11] F. Luca, P. Stanica, On the prime power factorization of n\, J. Number Theory 102 (2003), 
no. 2, 298-305. 

[12] B. Martin, C. Mauduit, and J. Rivat, Sur les chiffres des nombres premiers, submitted. 
[13] C. Mauduit, J. Rivat, La somme des chiffres des carrcs, Acta Math. 203 (2009), 107-148. 
[14] J. W. Sander, On the parity of exponents in the prime factorization of factorials, J. Number 

Theory 90 (2001), no. 2, 316-328. 
[15] T. Stoll, The sum of digits of polynomial values in arithmetic progressions, Functiones et 

Approximatio Commentarii Mathematici, accepted (Oct 21, 2011); preprint available from the 

author's webpagc. 

[16] G. Tenenbaum, Introduction a la theorie analytique et probabiliste des nombres, Belin, Col- 
lection Echelles, 2008. 

[17] W.-G. Zhai, On the prime power factorization of nl, J. Number Theory 129 (2009), no. 8, 
1820-1836. 

INSTITUT FUR DlSKRETE MATHEMATIK UND GeOMETRIE, TeCHNISCHE UnIVERSITAT WlEN, 
WlEDNER HAUPTSTRASSE 8-10, A- 1040 WlEN, AUSTRIA, 
E-mail address: johannes .morgenbesserOtuwien. ac . at 

Institut de Mathematiques de Luminy, Universite d'Aix-Marseille, 13288 Marseille 
Cedex 9, France, 

E-mail address: stollOiml.univ-mrs.fr 



